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Abstract 

The correlation energy in density functional theory can be expressed exactly in terms of the change 
in the probability of finding two electrons at a given distance rvx (intracule density) when the 
electron-electron interaction is multiplied by a real parameter A varying between (Kohn-Sham 
system) and 1 (physical system). In this process, usually called adiabatic connection, the one- 
electron density is (ideally) kept fixed by a suitable local one-body potential. While an accurate 
intracule density of the physical system can only be obtained from expensive wavefunction-based 
calculations, being able to construct good models starting from Kohn-Sham ingredients would 
highly improve the accuracy of density functional calculations. To this purpose, we investigate 
the intracule density in the A — > oo limit of the adiabatic connection. This strong-interaction 
limit of density functional theory turns out to be, like the opposite non-interacting Kohn-Sham 
limit, mathematically simple and can be entirely constructed from the knowledge of the one- 
electron density. We develop here the theoretical framework and, using accurate correlated one- 
electron densities, we calculate the intracule densities in the strong interaction limit for few atoms. 
Comparison of our results with the corresponding Kohn-Sham and physical quantities provides 
useful hints for building approximate intracule densities along the adiabatic connection of density 
functional theory. 
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I. INTRODUCTION 



Kohn-Sham (KS) density functional theory (DFT) (see, e.g., |lj) is a successful method 
for electronic structure calculations, thanks to its unique combination of low computational 
cost and reasonable accuracy. In the Kohn-Sham formalism, the total energy of a many- 
electron system in the external potential V ne = ^jfne(*"i) is rewritten as a functional of the 
one-electron density p(r), 



E[p) = T s [p) + U[p] + E xc [p] + J dvv ne (r) p(r) 



1.1' 



In Eq. (11. II) . T s [p] is the kinetic energy of a non-interacting system of fermions (usually called 
KS system) having the same one-electron density p of the physical, interacting, system. The 
Hartree energy U[p] is the classical repulsion energy, U[p] = ~ f dr f dr'p(r)p(r')\r — r'| _1 , 
and the exchange-correlation functional E xc [p] must be approximated. 

Despite its success in scientific areas ranging from material science to biology, DFT is 
far from being perfect, and a huge effort is put nowadays in trying to improve the approx- 
imations for E xc [p] (for recent reviews see, e.g., The focus of a large part of the 
scientific community working in this area has shifted from seeking explicit functionals of 
the density (like the generalized gradient approximations - GGA), to implicit functionals, 
typically using the exact exchange E x [p], which is only explicitly known in terms of the 
Kohn-Sham orbitals 0i(r) (for a recent review, see {4]). In this framework, although DFT 
was originally formulated as a method "without wavefunction" , it might be natural to go 
back and actually think of DFT approximations in terms of model wavef unctions. This way 
of thinking can be very helpful for building^ approximations, and for combining DFT with 
other many-body methods (see, e.g. 
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The adiabatic connection formalism (for a review, see 15[) is a useful tool to thin 
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DFT functionals in terms of wavef unctions. In its simpler and original version 
the electron-electron repulsion operator V ee in the iV-electron hamiltonian H (in Hartree 
atomic units used throughout), 

JV r,o N N 



H = T + V ee + V n 
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i=l 
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i>j=l 



1 



i=l 



[1.2) 



is multiplied by a real parameter A, which varies between and 1. At the same time, the 
external potential v ne (r) is replaced by another local potential, t> A (r), determined by the 



condition (allowed by the Hohenberg-Kohn theorems [19j], if p is v -represent able for all A) 
that the one-electron density p(r) do not change with A. In this way, we define a set of 
hamiltonians H x , 

H x = f + XV ee + V\ p A (r) =p(r) VA (1.3) 

all having the same p(r) as the one of the physical hamiltonian of Eq. (11.21) . In particular, 
at A = we have the KS hamiltonian, i.e., the hamiltonian of a non- interacting system 



of fermions with the same density of the physical system, and v 



A=0/ 



UKs( r )> the fa- 



miliar Kohn-Sham potential. If we denote by \I/ A the ground-state wavefunctions of each 
hamiltonian H A of Eq. (II. 3p , we easily find 



T s [p] 
EM 

EM 



(^ A=0 |T|^ A=U ) 



A=0\ 



{^ ={i \V ee \^ ={) )-U[p\ 



(^\V ee \^)d\-EM\-U[p}. / 



;i.4) 
;i.5) 
;i.6) 



where \1> A=0 is, in most cases, a single Slater determinant formed by the KS orbitals 0j. 

Equations (jl.5p and (11.61) can be rewritten in terms of the intracule density /(r 12 ) (also 
called in the DFT community spherically and system-averaged pair density), which was 
first introduced in the historical paper of Coulson and Neilson 20|]. Since then, /(?"i 2 ) has 
been used by several authors to understand electronic correlation both in density functional 



theory (see, e 
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22j|) and in post-Hartree-Fock methods (see, e.g., 



23 



24 



25 



30j). Given an iV-electron wavefunction the intracule density /(?"i 2 ) is defined 



as the integral of |\1/| 2 over all variables but r 12 = |r x — r 2 | 
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|^(ri2,R, r 3 , ...,rjv) 
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dHdr^...dr 
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where r 12 = r x — r 2 , and R = ^{ r i + r 2)- Here we have normalized /(ri 2 ) to the number 
of electron pairs. The quantity /(r 12 )47rr^ 2 is proportional to the probability distribution 
for the electron-electron distance in the state described by the wavefunction ^. Gill and 
coworkers 
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341 ] have defined an interesting "family of intracules", and made the 



hypothesis that the correlation energy of Hartree-Fock theory can be approximated as a 
linear functional of one of these intracules. 

In terms of the intracule densities J A (ri 2 ) associated to each wavefunction \I/ A of the 
adiabatic connection of Eq. (I1.3p . the Kohn-Sham correlation energy E c [p] of Eq. (11.61) can 



be rewritten exactly as 

'' I X (r 12 ) - I x=0 i 



Ec[p] 



o 



dX / dr 



ri2 
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rfA / rfri 2 47rri 2 [/ A (ri 2 ) - / 
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Correlation in Kohn-Sham DFT is thus fully determined by the change in the intracule 
density when the electron-electron interaction is turned on with the one-electron density 
p(r) fixed. The difference I A=1 (ri 2 ) — / A=0 (ri 2 ) determines the correction due to correlation 
to the expectation of V ee , and the integration over A recovers the correction to the expectation 
of T. By construction, there is no correction to the expectation of V ne . 

Starting from the observation that /(r 12 ) couples to the operator V ee in the same way as 
p(r) couples to V ne , i.e., that the expectations (\&|V^e|\I/) and (^/|V^=e|^) are given by linear 
functional of p(r) and 1(7*12), respectively, 



(*|Ke|*) = / dw ne {r)p{v), (Wee\*) 



dr 12 — I(r 12 ), 
r\2 



;i.9) 



it is possible to derive an exact formalism |35l . 1361 . 1371 . 1381 . 13 Ql ] in which a set of effective 
equations for each 7 A (ri 2 ) along the DFT adiabatic connection is coupled to the KS equations 
to generate the correlation energy from Eq. (11.81) . In this computational scheme one needs 
to make two approximations: 

1. the exact equation for I A (r 12 ) involves the solution of a many-body problem for a 



cluster of interacting fermions 



36], which is approximated with a radial Schrodinger 



equation for \pMjv2) 39J, possibly divided into effective geminals g A (r 12 ) 
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4Q|, 



N„ 



[-V* 12 +u&(n a )] 9i(ri2) = e A ^ A (r 12 ), J A (r 12 ) = \9i(rn)\ 2 , (1.10) 

i=l 

for which one needs to choose the number N g and the occupancy z^; 
2. an approximation for w A ff (ri 2 ) needs to be designed. 



As far as point 1 is concerned, we can say that the choice N g = 1 is always possible [39[, 
and yields good results in the uniform electron gas when combined with an approximation 
for w^ s (r 12) inspired to the Fermi- hypernetted-chain approach 4l|. Again in the case of the 
uniform electron gas, the choice of a determinant-like occupancy for the effective geminals 
(N g = N(N — l)/2, Vi = 1 (3) for even (odd) relative angular momentum states) yields 
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accurate results with much simpler approximations for w A ff (ri 2 ) [40|, |42j. In general, the 
choice of using localized geminals would make it easier to impose size consistency. 
Regarding point 2, the basic idea is to write u>e ff (V 12 ) as 



^eff( r 12) 



A 



Wes°(r 12 ) + — + w x c (r 12 ) 

fl2 
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The interaction Wgff°(r 12 ) is the one that, when inserted into Eqs. (11.101) . yields the intracule 
density of the KS system, / A=0 (ri 2 ), which can be constructed by inserting the KS Slater 
cl J1.7P - In this step, the analytical integrals developed by Gill and 



determinant into 
coworkers 
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34( to calculate Hartree-Fock intracules may turn extremely useful. 



In Eq. ( II. lip the term A/ri 2 ensures that the corresponding J A (ri 2 ) satisfies the electron- 
electron cusp condition (see, e.g., 43]). We need then to approximate w7 A (ri 2 ), an effective 
potential that should essentially "tell" to the intracule density that, while the electron- 
electron interaction is turned on, the one-electron density p(r) does not change. As the 
information on p(r) has been "washed away" in the integration over the center of mass R 
of Eq. (I1.7p . this constraint can be imposed only in an approximate way. For two-electron 
systems, for which Eq. (11.101) is exact with one geminal, simple approximations (based on 



the same ideas used in the uniform electron gas) for w A (ri 2 ) give accurate results |35l. 1371.138 1. 

To go one step further, that is being able to construct approximations that work for 
many-electron systems of nonuniform density, a crucial issue is to investigate the effect on 
I A (r 12 ) of the constraint of fixed p(r) as A increases. To this purpose, in this paper we 
address the following question: what happens to J A (r 12 ) when A — > 00 ? Although at first 
glance this question may seem purely academic, there are several reasons for investigating 
this strong-interaction limit of DFT. The intracule density of the physical system (A = 1) 
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can be obtained only from expensive wavefunction-based calculations (see, e.g 
and references therein), while in the A — > 00 limit the many-electron problem becomes 
mathematically simple, and we have recently shown 47j] that a solution can be constructed 
starting from the density p(r) only. The A — > 00 limit tells us what is the maximum extent 
to which the electrons can avoid each other without breaking the constraint of being in the 
given density p(r). This information can be very useful for constructing approximations. 
Last but not least, the strong interaction limit can be used to build interpolations between 
the non-interacting KS limit (A = 0) and the A — > 00 limit, yielding an approximation for 
the physical (A = 1) system [481 ] . 



This paper is organized as follows. In Sec. [TT]we derive and discuss the equations needed to 
calculate the intracule density in the strong-interaction limit (A — > oo) of DFT. In Sec. lIIIt we 
apply the equations of Sec. [IT] to calculate the A — > oo intracule densities of small atoms, by 
using accurate correlated one-electron densities p(r) as input. The results are then analyzed 
and discussed in Sec. HV] and the last Sec. |V]is devoted to conclusions and perspectives. 



II. THEORY 



The stron g-in teraction limit of DFT is defined by the A — > oo limit of the hamiltonians of 



Eq. O 
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501 ] . The mathematical details of this limit can be found in 47j. Here, we 



briefly summarize the physical ideas that lie behind the theory, only reporting the equations 
that will be used in the following sections. 

As A grows, the electrons repel each other more and more strongly. However, they are 



forced by the externa 
can be shown 
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potential V x of Eq. (11.31) to yield the density p(r). As A — ► oo, it 



50 ] that, in order to keep the electrons in the density p(r), V x must 
be proportional to A, V x -*°° — » \V, In this limit, the kinetic energy becomes negligible (of 
orders \/A \^), and the solution of H x ^°° reduces to a classical equilibrium problem for 



the 3N dimensional function 



N 



E 



pot 



>!,..., r N ) = 



N 



V(Ti 



v(r) = lim 



v x (r) 



A— >oo A 



(2.12) 



i>j=l 1 J i=l 

The square of the corresponding wavefunction, |v|/ A ^ 00 | 2 ) becomes a distribution that is zero 
everywhere except in the configurations (r^\ ...r^) for which £ , pot (r 1 , r N ) has its absolute 
minimum. Typically, for a reasonable attractive potential v(r), E pot has a discrete set of 
minimizing configurations. In this case, however, the density corresponding to |\1/ A ^ 00 | 2 
would be given by a sum of peaks centered in the minimizing positions rf \ p(r) oc Yli — 
). In order to get a smooth density like the one we find in the quantum mechanical 



problem at A = 1, we need a special potential v(r) in Eq. (12.121) : the potential v(r) must 
make the minimum of the 3A-dimensional function E pot degenerate over the 3-dimensional 
subspace M j^] 

M = {r 1 = r, r 2 = f 2 (r), r N = f N (r), r G P}, (2.13) 

where P is the region of space in which p(r) ^ 0. From the physical point of view, the 
distribution |\I' A ^ 00 | 2 ) which is zero everywhere except on M, describes a state in which the 



position of one of the electrons can be freely chosen in P, but it then fixes the positions of 



3- 



all the other N — 1 electrons via the co-motion functions f,;(r) |47j. The strong- interaction 
limit of DFT is thus the generalization of the more familiar Wigner-crystal state to smooth 
densities. In the Wigner crystal state, in fact, the constraint of having a given density is 
relaxed, and p(r) becomes typically proportional to £\ 5(r — r^), losing any resemblence 
with the quantum mechanical A = 1 density of atoms and molecules. 

From the condition that E pot have its minimum over the entire subspace M and that 
the electrons be indistinguishable, one finds that the co-motion functions fj(r) must satisfy 
special properties, which are reported in [42 1. To determine the co- motion functions from 
the density p(r), we use the observation |47| that, since the position of the first electron 
determines the positions of all the others, the probability of finding the first electron in the 
volume element dr around the position r must be the same of finding the i th electron in the 
volume element c?fj(r) around the position fj(r). This means that all the co- motion functions 
fj(r) must satisfy the differential equation 

p($(r))dfi(r) =p(r)dr, i = 2, ...,JV. (2.14) 

In order to construct the co-motion functions we thus have to find the initial conditions for 



the integration of (12.141) which also satisfy the properties reported in 471 ] . An example of 
such calculations for spherical densities is carried out in j^j]. The strong interaction limit of 
DFT is thus entirely determined by the co- motion functions fj(r), which can be constructed 
from the density via Eqs. (12.141) . 

To obtain the intracule density I x ~^°°(ri2) corresponding to the distribution jil/-^ 00 ! 2 , we 
have to consider that the electron-electron distance only depends on the position of the first 
electron, r. By defining the N(N — l)/2 distances dij(r) for which |v|/ A ^ 00 | 2 is non zero, 

rf i ,(r) = |f,,(r)-f i (r)|, i,j = l,...,N, K (2.15) 

I — I 

(with fi(r) = r), and by considering that each position r has a probability ^p(r) [47j, we 
obtain 

4 7 rr? 2 P Woo (n 2 ) = J2 fdv^-Sin.-d^v)). (2.16) 

i>j=i J 
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III. APPLICATION TO ATOMS 



We consider here the case of spherical densities, and we apply Eq. (12.161) to few atoms. 
When pjr) = p(r), the A — > oo problem can be separated into an angular part and a radial 
part [471 ] . The distance r from the nucleus of one of the electrons can be freely chosen, and 
it then determines the distances from the nucleus of all the other N — 1 electrons via radial 
co-motion functions /i(r), as well as all the relative angles afy-(r) between the electrons j^ . 
The radial co- motion functions fi(r) can be constructed as follows [47]. Define an integer 
index k running for odd N from 1 to (N — l)/2, and for even N from 1 to (N — 2)/2. Then 



Mr) 



N~ l (2k- N e (r)) r < a 2k 
N~\N e {r)-2k) r>a 2k 



, N-\N e (r) + 2k) r<a N _ 2k 

f2k+i(r) = < , (3.17) 



N-\2N -2k- N e (r)) r > a N _ 



2k, 



where Oj = N e 

N e {r)= I A7Tx 2 p(x)dx, (3.18) 
Jo 

and N~ 1 {y) is the inverse function of N e (r). For odd N, these equations give all the needed 
N — 1 radial co-motion functions, while for even iV we have to add the last function, 

f N (r) = N-\N-N e (r)). (3.19) 

The relative angles ctij(r) between the electrons can be found by minimizing numerically the 
electron-electron repulsion energy 2~2i>j(fi( r ) 2 + fj( r ) 2 ~ ^fi( r )fj( r ) cosa^-) -1 / 2 . The radial 
co-motion functions of Eqs. (13 .171) - (13. 191) satisfy Eq. (I2.14|) for spherically symmetric p, 

4tt /,(r) 2 p(/ i (r)) \f[{r)\ dr = 4vr r 2 p(r) dr, (3.20) 

and, together with the minimizing angles a>ij(r), yield the minimum of E pot of Eq. (12 . 1 2[) for 
spherically symmetric v(r) [47j]. Physically, the solution of Eqs. (13. 171) - (13. 191) makes the N 
electrons always be in iV different spherical shells, each of which contains, on average in the 
quantum mechanical problem (at A = 1), one electron. In the A —>■ oo limit, the electrons 



become strictly correlated, and all fluctuations are suppressed (see, e.g., [51]): the space is 



divided into iV regions, each of which always contains exactly one electron. 
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0.5 1 1.5 2 2.5 3 3.5 4 
r 

FIG. 1: The electron-electron distance d±2 for the He atom density in the strong-interaction limit 
of DFT. d\2 is, in this limit, completely determined by the distance r from the nucleus of one of the 
two electrons, d±2 = di2(r). The value r = a±, for which d\2{r) has its minimum, corresponds to 
the radius of the sphere containing, on average in the quantum mechanical problem, one electron, 
Jq 1 4irr 2 p(r) dr = 1. All quantities in Hartree atomic units. 



For spherically symmetric densities, the electron-electron distances of Eq. (12.151) then 
become 

dy(r) = ^(r) 2 + f,{rf - 2/ i (r)/ j (r) cosc^(r). (3.21) 

We can invert Eq. (I3.2ip (if dij(r) is non monotonous we have to invert each monotonous 
branch separately), and simplify Eq. (I2.16P into 

4vrr 2 2 /^(r 12 ) = £ p(d^(r 12 )) |^-(^( ri2 ))|-i, p ( r ) = (3 . 22) 
i>j=i 



As a starting point, it is instructive to analyze the simple case of the He atom. Here, 
we have used the accurate variational wavefunction of Ref. (52] (see also 53] and 35]) to 
generate the density p(r), from which / A ^°°(ri2) is constructed. From the same accurate 
variational wavefunction, we have also computed the intracules of the KS system (see 35]), 
/ A=0 (r 12 ), and of the physical system, J A=1 (r 12 ), which will be compared to / A ^°°(r 12 ). When 
N = 2, in the A — ► 00 limit the relative angle between the electrons becomes always a%2 = 7r 
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FIG. 2: The intracule density of the He atom along the linear adiabatic connection of DFT: the 
three intracules correspond to three systems with the same one-electron density p(r) and electron- 
electron interaction \jr\ 2 . The intracules at A = (KS system) and A = 1 (physical system) have 



53| and 



been constructed from the accurate variational wavefunction of Ref. [52| (see also 
The intracule at A = oo is obtained in this work, as described in the text. All quantities in Hartree 
atomic units. 



(maximum angular correlation), and we have only one co-motion function, f 2 {r), which fully 



determines the electron-electron distance 



47 



5fJ, 



f 2 (r) = N-\2 - N e (r)), d 12 (r) = r + f 2 (r). 



(3.23) 



The function di 2 (r) is reported in Fig. [TJ It has a minimum for r = a± — N~ l (l), the radius 
of the shell containing, on average in the quantum mechanical problem, one electron. We 
have f 2 {a\) = cq, and di 2 (ai) = 2cq. The two electrons, thus, never get closer than 2cq, so 
that f A ^°°c ^ 2a,i) = 0. From the indistinguishability of the two electrons, we have the 
property 



47 



501 ] f 2 (f 2 (r)) = r, which, when combined with Eq. (13.201) . shows that the two 



invertible branches of <i 12 (r) (corresponding to < r < cq and to r > cq) give the same 
contribution to Eq. ( I3.22H . We can thus just invert the function di 2 (r) in r G [0,cq] and 
multiply the result by 2. This is a general property, also valid for N > 2 [47]: we can always 
invert the functions dij(r) in the domain r £ [0,cq], and then multiply the result by N. In 
Fig. [2] we report the intracule densities of the He atom for A = (KS), A = 1 (physical) and 
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FIG. 3: The electron-electron distances dij for the Li atom density in the strong-interaction limit 
of DFT. The distances dij are, in this limit, completely determined by the distance r from the 
nucleus of one of the electrons, dij = dij(r). The value of r is varied here between and cq, 
the radius of the sphere containing, on average in the quantum mechanical problem, one electron, 
Jq 1 Anr 2 p(r) dr = 1. As explained in the text, the case r > a\ does not need to be considered, 
since it simply corresponds to interchanging two or more electrons, yielding the same values for 
the electron-electron distances. All quantities in Hartree atomic units. 



A — > oo. Because d' 12 (ai) = 0, / A ^°°(r 12 ) has an integrable divergence when r 12 — > 2a\ 

K«i) S/a 



p(r) = 2nr p(r). 



(3.24) 



n2->2a+ v-j/(ai)(ria - 2ai) 
Divergences come from the fact that the quantities we calculate here are the distributions 
towards which the physical quantities tend when A — > oo. This aspect is clarified with a 
simple example in the Appendix of Ref. {47]. 

Using Eqs. fl3.17l) - (l3.19p . we have calculated the electron-electron distances dij{r) of 
Eq. (13.211) for Li, Be and Ne. For the Li atom, we have used the fully correlated den- 
sity of Bunge 54|, and for the Be and the Ne atoms the accurate densities of Ref. 55]. All 
calculations are done numerically, on a grid. In Figs. [3] and HJ we report dij(r) for Li and Be 
when < r < a\. Electrons are labeled with numbers 1, 2, 3,..., meaning that electron 1 is in 
the shell < r < a\, electron 2 is in the shell cq < f^fr) < 0-2, and so on. Correspondingly, 
the distances are labeled 1-2, 1-3, etc. As discussed for the case of the He atom (for further 
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FIG. 4: Same as Fig. [3] for the Be atom density. 
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details see 47j|). we only need to consider N times this situation, since exchanging two or 
more electrons always correspond to the same physics, and thus to the same values of the 
electron-electron distances. 

In Fig. [5] we show the intracule density multiplied by the volume element 4:irrf 2 at A = 0, 
A = 1 and A — > oo for the Be atom. The case of the Ne atom is displayed in Fig. El where 
we show the intracule densities multiplied by Anr^: the area under each curve gives the 
expectation (\l/ A |V r ee |\I/ A ). In both cases, the KS intracules at A = have been constructed 



from the accurate Kohn-Sham potentials of Ref. [55j, while the intracules of the physical 
system (A = 1) are obtained from variational Quantum Monte Carlo results [44J. Integrable 
divergences in the A — > oo intracules appear, as for the He atom case, at elecron-electron 
distances for which c^(r) = 0. Notice that, due to the strict correlation at A — > oo, there 
is a finite minimum distance r™ in > between any pair of electrons, so that I x ~^°°(r 12) = 
when r\2 is less than the smallest of the r™ m . 



IV. DISCUSSION OF RESULTS 



The A — > 00 limit of DFT describes the case of maximal angular correlation and maximal 
radial correlation between the electrons, compatible with the constraint that the probability 
of finding one electron at postion r be equal to p(r)dr, where the density p(r) corresponds to 
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FIG. 5: The intracule density of the Be atom multiplied by the volume element 47rrf 2 along the 
linear adiabatic connection of DFT: the three intracules correspond to three systems with the 
same one-electron density p(r) and electron-electron interaction A/ri2- The intracule at A = (KS 

55[ and the intracule at A = 1 



system) has been obtained from the accurate KS potential of Ref. 
(physical system) is taken from Ref. [44|. The intracule at A = oo is calculated in this work, as 
described in the text. All quantities in Hartree atomic units. 



the quantum mechanical hamiltonian of Eq. (jl.2p . For atoms, we see from Figs. [21 [5] and [6] 
that the intracule of the physical system (A = 1) is much closer to the KS intracule (A = 0) 
than to the A — > oo limit, as expected for weakly correlated systems. In more correlated 
situations like streched bonds, we can expect the A = 1 case to be more in between the 
A = and the A — > oo limits. The investigation of such cases will be the object of future 
work. 

What can we learn from the intracule densities in the A — > oo limit? As a first step, we 



have reconsidered the results of Ref. [35j for the He-like ions. In that work, we had shown 
that a good approximation for of Eq. (11. lip was given by the screening potential of a 
sphere of uniform density, charge 1, and radius r s to be determined, 

w^r 12 ) « - 1 / dr. (4.25) 

\ / J\r\<r a l r ~ r 12| 



This approximation was first introduced by Overhauser for the uniform electron gas 56], 



where r s was set equal to the usual density parameter r s , i.e., the radius of the sphere con- 
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FIG. 6: The intracule density of the Ne atom multiplied by Airr\2 along the linear adiabatic 
connection of DFT: the three intracules correspond to three systems with the same one-electron 
density p(r) and electron-electron interaction X/ri%. The intracule at A = (KS system) has been 



obtained from the accurate KS potential of Ref. 
is taken from Ref. 



551 ] and the intracule at A = 1 (physical system) 



441 ] . The intracule at A = oo is calculated in this work, as described in the text. 
The area under each curve gives the expectation (\E ,A |V^e|\E' A )- All quantities in Hartree atomic 
units. 



taining, on average, one electron. In Ref. [40] it was shown, by comparison with Quantum 
Monte Carlo results, that for the uniform electron gas the Overhauser potential very accu- 
rately recovers the short-range part (r\2 < r s ) of /(r^). For the He-like ions, setting 35] 



where p is an average density, P = f rfrp(r) 2 , also yields accurate results for 



the short-range part of /(r 12 ), with "on-top" values 1(0) essentially indistinguishable from 
those coming from very accurate Hylleras-type variational wavef unctions. Here, we realized 



that the values r, used in Ref. 



351 ] for the He-like ions are, within few percents, equal to 



the values cq of the radii of the sphere (centered in the nucleus) containing on average one 
electron, thus making Eq. (I4.25[) work equally well for the uniform electron gas and for the 
He series, with the same choice for the screening length r s . As shown in the previous sec- 
tion, the value cq plays a special role in the A — > oo limit of the He-like ions. We can thus 
hope to learn from the A — > oo intracules something about "multiple screening lenghts" for 
approximating (r^) for many-electron systems of nonuniform density. 
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FIG. 7: Be atom: the core-core, core- valence and valence- valence contributions to the intracule 
densities at A = and A = oo of Fig. [5j 
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FIG. 8: Ne atom: the core-core, core- valence and valence- valence contributions to the intracule 
densities at A = and A = oo of Fig. [6j 
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To this purpose, and with the idea in mind that to have a size-consistent method we need 
to use localized geminals in Eqs. fll.lOp . we have further analyzed our results by dividing them 
into core-core, core-valence, and valence-valence contributions, comparing the A — > oo and 
the A = case. We consider here the Be and the Ne atoms. For the A — > oo case, the core- 
core contribution comes from the distance 1-2, corresponding to the two electrons that are 
in the sphere containing, on average in the quantum mechanical problem, 2 electrons. The 
distances of electrons 1 and 2 from the other electrons define the core-valence contribution, 
and the rest is the valence- valence part. For the KS system (A = 0), we have orbitals, so the 
three contributions are defined in the usual way, using the quantum mechanical shells (Is 2 
for core-core, etc.). The three contibutions are shown in Fig. [7] for Be and in Fig. [S]for Ne. 
In the core- valence case of Be, we see that the extremely correlated A — > oo limit differs from 
its KS counterpart only in the short-range part, r\2 < 2.5. The valence-valence case of Be 
resembles the two-electron case of the He atom. The core-valence and the valence-valence 
contributions to the Ne atom also show, essentially, correlation of short-range type, even in 
the extreme A — > oo case. In our future work, we plan to use these results to build and test 
approximations for w^{r\2}. 



V. CONCLUSIONS AND PERSPECTIVES 



We have calculated, for the first time, the intracule densities for small atoms in the strong- 
interaction limit of density functional theory. Our results can be useful to better understand 
correlation in the density functional theory framework, and to build approximations for 
correlation energy functionals based on intracules. Our future work on this subject will 
address several points: 

1. The generalization of this calculation to non-spherical densities; 

2. The study of the next leading term in the A — > oo limit, thus including zero-point- 
motion oscillations; 



3. The use of these results to fully develop the ideas of Refs. [35j, |36|, |37|, |38|, |39|, in which 
an approximation for the correlation energy in density functional theory is contructed 
from effective equations for the intracule density. 
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